arXiv:hep-ph/0112193vl 14 Dec 2001 


Electron-positron pair production in the external electromagnetic 
field of colliding relativistic heavy ions 


A. Aste^, G. Baur^, K. Hencken^, D. Trautmann^ and G. Scharf^ ^ 

^ Institute for Theoretical Physics, University of Basel, Klingelbergstrasse 82, 4054 Basel, Switzerland 
^ Institut fiir Kernphysik (Theorie), Forschungszentrum Jiilich, 52425 Jiilich, Germany 
® Institute for Theoretical Physics, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland 

November 30, 2001 

Abstract. The results concerning the e''"e“ production in peripheral highly relativistic heavy-ion collisions 
presented in a recent paper by Baltz et al. are rederived in a very straightforward manner. It is shown that 
the solution of the Dirac equation directly leads to the multiplicity, i.e. to the total number of electron- 
positron pairs produced by the electromagnetic field of the ions, whereas the calculation of the single pair 
production probability is much more involved. A critical observation concerns the unsolved problem of 
seemingly absent Coulomb corrections (Bethe-Maximon corrections) in pair production cross sections. It 
is shown that neither the inclusion of the vacuum-vacuum amplitude nor the correct interpretation of the 
solution of the Dirac equation concerning the pair multiplicity is able the explain (from a fundamental 
point of view) the absence of Coulomb corrections. Therefore the contradiction has to be accounted to the 
treatment of the high energy limit. 
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1 Introduction 

The problem of e“'"e“ pair production by the collision of 
highly relativistic nuclei has attracted a lot of interest dur¬ 
ing the past few years, since this process will have implica¬ 
tions for future experiments performed at new facilities as 
BNL’s RHIG and CERN’s LHC. From the theoretical side, 
the result that the Dirac equation can be solved exactly 
in the electromagnetic background field created by two 
nuclei in the limit where the two nuclei are ’ultrarelativis- 
tic’ [Ql^, seemed to lead to the unexpected consequence 
that the single pair cross section is equal to its Born value 
(given by the second order diagram in Fig. ^) [|. The term 
’ultrarelativistic’ is used here for the limiting case where 
the Lorentz factor 7 becomes large (7 —> 00), i. e. the ve¬ 
locity of the colliding ions approaches the speed of light. 
Of course the problem of pair production in an external 
held goes back to the beginning of QED [^,|^. In connec¬ 
tion with the relativistich heavy ion colliders it was found 
that the impact parameter dependent probability in per- 
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^ Note that this is only true for the total cross section. As 
it was shown in 0, this is not true for impact parameter de¬ 
pendent probabilities. 


turbation theory can become larger than one, which was 
shown to result in multiple pair production [ 

A series of papers on the same topic followed 00: 
with the aim to show that there must be an error 
in the interpretation of the results found in ||^. An inter¬ 
esting observation was hnally made in a paper by Baltz et 
al. p^ , where it was shown that the expression for pair 
production derived from the results in |^] describes the 
total number of produced pairs, and not the single pair 
production. Additionally, the importance of taking into 
account the vacuum-vacuum transition amplitude when 
going over from wave mechanics to the full external held 
problem with quantized electron held was pointed out, as 
was already discussed in |||,^ . A recent work which treats 
the structure of Goulomb and unitarity corrections to the 
single and multiple pair production is |p^ . 

We will derive in this paper the main results obtained 
in fl^ in a very compact way using a different held theo¬ 
retical point of view. The correct expression for the mul¬ 
tiplicity can be derived indeed in a very straightforward 
manner from the fundamental equations which dehne the 
S-matrix for the external held problem. We establish the 
connection between the single-particle matrix elements 
from the solution of the Dirac equation and the pair pro- 
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duction in a full many-body theory. Furthermore, we give 
a strong argument which highlights that the results pro¬ 
vided from the theory of fermions in a classical electro¬ 
magnetic background field are insufficient to explain the 
dicrepancies concerning Coulomb corrections in the liter¬ 
ature. 

It is therefore mandatory to understand the correct 
interpretation of the results obtained from ’infinite 7 ’ cal¬ 
culations and their implications on real electromagnetic 
processes, as it was suggested in the work of Lee and Mil- 
stein [|^ where they have shown that Coulomb corrections 
can be obtained by the correct regularization of matrix el¬ 
ements. 


2 The external field approximation for QED 

For highly relativistic heavy-ion collisions one usually as¬ 
sumes that the ions of electric charge Zie and ^26 are suf¬ 
ficiently energetic and massive so that the deviation from 
straight-line trajectories can be neglected. The trajecto¬ 
ries of the two ions are then, respectively, z = ±t {c = 1 ), 
xj^ = b/ 2 , and there is a gauge in which the potential of 
the two ions has the following form ||^ 

A^(t,x) = Zied{v+x)v^ ) 

-f log(^^—, ( 1 ) 

where = (1,0, 0, ±1)/v^. The fact that the two electro¬ 
magnetic fields are compressed due to the Lorentz contrac¬ 
tion to a sheet in the b plane, together with their structure 
in space-time allows for the exact solution of the Dirac 
equation in this case. 

We will deal in the following with the general external 
field problem without referring to the special form of the 
electromagnetic field created by ultrarelativistic ions as 
given by Eq. (|^), in order to clarify how the solution of 
the Dirac equation for the electron wave function and the 
full S-matrix for the quantized electron field in an external 
electromagnetic field are related to each other. 


and define a unitary S-matrix S = Since we pre¬ 

suppose free dynamics for t ± 00 , we can base second 
quantization on the Fock representation of the free Dirac 
field. The one-particle sector of the Fock space then con¬ 
sists naturally of a positive spectral subspace (electrons) 
and negative spectral subspace (positrons) of the Hilbert 
space of Dirac four-spinors 

Hi = (L2(7^3))^ = © P^Ui. (5) 

with the standard inner product. Here, denote the cor¬ 
responding spectral projection operators of the free Dirac 
Hamiltonian. The quantized Dirac field is then defined as 
usual in the Schrodinger picture 

V’(/) = KPlf) + d(P°/)+ = b[U) + d(/_)+. (6) 

The electron annihilation operator b{fj^) destroys an elec¬ 
tron with wave function /_|_, and the positron emission op¬ 
erator fi(/_)+ produces a positron with wave function /_. 
The annihilation and emission operators fulfil the fermionic 
anticommutation relations 

{b{f+),H9+)~^} = if+,9+), {d{f-y,d{g-)} = (/-,ff-), 

( 7 ) 

and the Fock vacuum 17 is defined by 

b{f+)n = d{f.)n = 0 V/ e Hi. (8) 

Note that &(/+) depends antilinearly on /+, whereas d{f-) 
depends linearly on /_. This explains the slight difference 
in the commutation relations for electron and positron 
operators in (|^). We mention that naturally, the widely 
used field operators can be defined which correspond to 
the operators introduced above according to 


Hur = j 

1 d^xb{x)+f+{x), 

(9) 

dif-r = 

[ d^xd(x)+/*(x) 

( 10 ) 


or, given an arbitrary orthonormal basis fj in P^TLi and 
Qk in P°Hi, 


2.1 Construction of the S-matrix 


In order to construct the full S-matrix of QED in the ex¬ 
ternal field approximation, we start from a one-particle 
dynamics, defined by a time-dependent Hamiltonian act¬ 
ing on a state described by a Dirac spinor 

H{t)=Ho + Vit), (2) 

where V(t) contains the interaction with the external elec¬ 
tromagnetic field given by 

H(t,x) = e((^(t, x) — aA(t, x)). (3) 


The potentials are assumed to vanish for t ±oo in such 
a way that the wave operators exist as strong limits 


Win = lim C/(t, 0 )+e-*-^“* 

ont t^±00 ' ^ 


='^digk)+gk{:>c) 

3 k 

( 11 ) 

such that = b{x.) + (i(x)+ fulfils the well-known dis¬ 
tributional identity 

{•0(x),?/'(x')+} = l4(5^^^(x-x'). (12) 

The one-particle S-matrix, obtained via its perturba¬ 
tion series 


5* = 


n =0 


+ 00 ti tn-l 

{-ir J dh J dt2... J dHH(ti)...H(t„) 


where 


( 4 ) 


V{t) = 


(13) 

(14) 
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or by solving directly the Dirac equation with the external 
field, can be written as 



if we define S±± = P^SP^ and = P^S+P^. It is 

often said that 5++ and S are related to electron and 

positron scattering, whereas _and S _^ are related to 

pair creation and annihilation. But the connection of S±± 
to different measurable quantities has to be treated care¬ 
fully. It should be stressed that the one-particle Dirac the¬ 
ory has no consistent physical interpretation, but only the 
corresponding second quantized many-particle theory in 
Fock space which we are now going to discuss. 

If one lifts the one-particle Hamiltonian with indefinite 
energy spectrum to the Fock space 

H =J2[{fj,Hf,)b{f,)+b{h) - {gj,Hgk)dig,)+dig ,)\, 

jk 

(16) 

then the Dirac field operator fulfils the Heisenberg equa¬ 
tions of motion 




(17) 


= [^{U ^(t,0)/)+,H]. 


(18) 


Roughly speaking, ( p^ ) simply accounts for the fact that 
the total energy of a general Fock state should be given 
by the sum of energies of each single particle. Note also 
the important minus sign in (16). 

The definition (^) of the Fock vacuum would lead in 
the case of time-dependent external fields to Fock repre¬ 
sentations which continuously change in time. To avoid 
this somewhat odd construction, we retreat to scattering 
theory by defining the second quantized S-matrix S in cor¬ 
respondence to the Heisenberg equations of motion ( p^ 

^{S+f) = S-V(/)S, (19) 

^(5+/)+ = S-V(/)+S V/eHi. (20) 

It is even possible to give an explicit normally ordered 
expression for S 0,0 


5 ^fj^s+-szlb+d+ . ^{s+-^-i)b+b . 


. g(l-‘STi)iid+ . + 


( 21 ) 


where the terms in the exponents are a shorthand for the 
operators 


Ab+b = Y.{f„Afk)b{fi)+b{h), (22) 

jk 

Ab+d+ = Y,ifj,^9k)bifj)+d{gk) (23) 

jk 

and similarly for Add^ and Adh. Note that Sz\ is the in¬ 
verse of the matrix S _alone and therefore not identical 


to = P^S~^P^. C denotes the vacuum-vacuum tran¬ 
sition amplitude (12, Sf7). This vacuum-vacuum transition 
amplitude is always present, but whereas it has absolute 
value one, e.g., for static fields (and is therefore dropped 
in the calculations), C is indispensable for the unitarity of 
S, as has been verified in [0 by perturbation theory. In 
our case where two static Coulomb fields are moving rel¬ 
atively to each other, the lowest order contribution to the 
vacuum-vacuum transition amplitude is given by a loop 
graph shown in FigJ^. It is a well-known result that all 
vacuum graphs appearing in pair production amplitudes 
(see the example in Fig. ||) can be absorbed in the factor 
C. Additionally, the lowest order contribution is finite due 
to gauge invariance, whereas higher order contributions to 
the vacuum-vacuum transition amplitude are finite due to 
power counting arguments. Therefore, the theory is free 
of any regularization problems. 



Z. 



Fig. 1. Lowest order contribution to the vacnum-vacuum tran¬ 
sition amplitude and single pair production graph containing 
a vacuum bubble. 


Equation is the important link showing the rela¬ 
tion between the (single-particle) S matrix from the so¬ 
lution of the Dirac equation and the general many-body 
theory. 

We finally mention that the S-matrix S in Fock space 

exists, if and only if _is a Hilbert-Schmidt operator. 

This condition is not satisfied with the potential (|^) as it 
stands, but it can be achieved by appropriate regulariza¬ 
tion. 


2.2 Multiplicity 

From (H) it follows immediately 

6 (/)S = S[6(5++/) + d(5l+/)+], 


(24) 
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a simple identity which allows to derive easily the average 
number of electron-positron pairs produced by the exter¬ 
nal field. The number of pairs in the final state is equal 
to the number of electrons in the final state, which can be 
calculated from the electron (or positron) number opera¬ 
tor 

7V, = ^&(/„)+6(/n), fVp = ^d(g„)+d(5„). (25) 

n n 

Therefore, we obtain for the number of pairs in the final 
state created by the external field out of the vacuum SJ7 

iV= (Sf2,^6(/„)+6(/„)Sf2) 

n 

= J2iHfn)sn,b{fn)sn). (26) 

n 

Using (^J), this can be written as 

N = '£{d{sl+fn)n,disl+u)n) = '£{fn,s+-sl+u) 

n n 


The matrix Ai = S'-]_ S_]_ contains the additional scat¬ 

tering factor SZ -, which is missing in the literature when 
differential cross sections for pair production are consid¬ 
ered. As will be explained below, Ai is identical to the 
result using the usual Feynman propagators. 

Unfortunately, the calculation of an integral operator 
like SZ- is not a trivial task, but from Eq. (H) we can 
derive the unitarity condition 

S-+S+_ + S—Sl:_ = P° , (33) 

which can be multiplied from the left by 5"^_ SZ- leading 

to the equation 

Ai = S+-S+_ + AiS-+S+_ , (34) 

which can be iterated. 


2.4 n-pair production 

The vacuum-vacuum transition probability jUp can be 
derived from 


= Tr(5+_5^+) (27) 

or equivalently, using the positron number operator 

N = Tr{S-+S+_). (28) 

Therefore, the results obtained from the solution of the 

Dirac equation for S _j. [Q must not be misinterpreted as 

the single pair production amplitude. Furthermore, this 
result for which a lengthy derivation has been given in 
ijl^, follows here in an almost trivial way. 


1 = (Sf?, Sf2) = (35) 

with the result Q 

ICr^ = det(H-AT+A4). (36) 

It is clear that the probability Ptot to produce an arbi¬ 
trary number of pairs is equal to 1 — |Cp. The well-known 
identity 

logdet(l-f A4+A4) = Trlog(l-)-A4+A4) (37) 


2.3 Single pair production probability 


In order to derive the production probability for a single 
electron-positron pair with wave functions <Pe and we 
have to compute first the S-matrix element 

Mep= (6(<I>e) + d(^>p) + f2,Sf2) 

= C{fI,d{<Pp)+b{<P^)+e^+-^--^*‘^^ fl) . (29) 

The definition ( p^ and the anticommutation rules (j^) en¬ 
able us after some algebra to derive the commutation rules 

0 

= e^+-^"'^‘'^(6(/) + d((5+_5li)+P°/)+) , (30) 

= e^+-^--'^^<^^idif)-b{S+.SzlP°f)+) , (31) 

which reduce the matrix element to the form 


Mep = C(I2,d(<?p)d((5+_5Zi) + P°^>e)+D) 
= C{<Pe,S+-SZl<Pp). 


leads to 


|C| 


1 ( T w 

= exp -Tr 2^ - 


n—1 


n 


or 


1 ^ „ 1 . / ^ ^ (-Ai+Aiy 

.2^P„ = _-l = exp(-Tr2^-Z- 


|C |2 ^ 

' ' n—1 


(38) 


- 1 . 


. . CTT . 

Utilizing the expansion in ( pQl) allows to write down the n- 
pair production probability by collecting all terms in 
containing n Ai^Ai terms, e.g. 


Pi = 


|C|^ 


[(Tr(A4+A4))2 - Tv{{Ai+Aif)]. (40) 


This can be understood if one takes a look at the actual 
calculation of the n-pair production probability analogous 
to the single pair production calculation presented in sect. 
2.3. It follows also from the expression of the determinant 
in terms of antisymmetric tensor products which can be 
found in lIT^ 


det(l + A) = ^Tr(A"(A)), 

n—Q 


(32) 


(41) 
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if A is a trace-class operator. A” (A) is the tensor product 
of A over the antisymmetric n-particle space. 

If ’exchange terms’ of the form Tr((AI''"AI)") are ne¬ 
glected in (|39[), then the probability of producing n pairs 
goes over into a Poisson distribution 


p . ,-Tr(At+At) mM+M)r 
n\ 


(42) 


which seems to be a reasonable approximation due to the 
correlation between the momenta of electron and positron 
in this case, as was also found in the calculations per¬ 
formed in . The poisson distribution had already been 
derived theoretically in an earlier work . There, the sud¬ 
den (or Glauber) approximation and a quasiboson approx¬ 
imation for e+e“ pairs was assumed, leading to a sim¬ 
ple but very instructive example which elucidates the pair 
production mechanism. It is remarkable that the Poisson 
distribution gives the multiplicity of particles to be 

Npo^sson = Tr(M+M) (43) 

quite similar to the exact one of Eq. (^). 


2.5 Perturbation theory 


We briefly mention here the perturbative expansion of the 
previously discussed matrix elements. As it follows from 
the solution of the Dirac equation the perturbative expan¬ 
sion of S±± is given by retarded propagators 




^+m 

p2 _ ^2 _|_ jpgO’ 


(44) 


S'i”2(p,q) = -» ^^^p-i -P±(p)7°/ d^Pi 

4ip-pi)S'^‘'\pi)i{pi-p2) ...S'^‘'\pn-l)4iPn-l-q)PU^)- 

(45) 

The projection operators are given in momentum space 

by 

■P±(P) = ^(-E'±ap + /3TO): (46) 

where a and j3 are standard Dirac matrices. The pair pro¬ 
duction amplitude, on the contrary, has to be calculated 
using Feynman propagators 


S^{p) 


4+m 

p2 _ _|_ jQ ’ 


(47) 


(5+_5li)(")(p,q) 

= -^ (2^)2n-l -P+(P)b° / Al ... An -1 

4(p-Pi)A(pi)4(pi -P 2 ) ...S^{pn-l)4iPn-l - q)P°iq). 

(48) 

Eqs. ( ^5|J^ ) represent well-known results from (quantum) 
field theory which are discussed in many standard text¬ 
books I ^ ^ . An explicit calculation which shows how the 


transition from retarded to Feynman propagators takes 

place when and {S^ _S'zi)^"’) are compared can be 

found in . It is the inverse of S _which is responsible 

for this transition. 

Since it is possible to solve the Dirac equation exactly 
in the ultrarelativistic limit, it is clear that the retarded 
propagator is also known in this case, and is therefore a 
much simpler object than the Feynman propagator from 
the calculational point of view. 


3 Coulomb corrections and multiplicity 

It was argued in Q that the absence of Coulomb correc¬ 
tions in the high energy limit as derived in [Q could be 
attributed to the difference in the interpretation between 
the multiplicity N = ^(5*^_A+) and the single pair pro¬ 

duction probability Pi = |C'pTr(A4+A4). Here we want 
to present some arguments that show that such an expla¬ 
nation cannot hold. Coulomb corrections and ’unitarity 
corrections’ (following the terminology used in H) can 
be treated as two different phenomena, which stem from 
different sources. Whereas unitarity corrections are due 
to many-particle aspects, Coulomb corrections affect the 
matrix element A4 directly. 

Since the theory of pair production developed above is 
valid quite generally, we can consider the special case with 
two ions that have charge numbers Zi and Z 2 , in the limit 
where Zia remains fixed but Z 2 CX goes to zero. Expanding 
in powers of Z 2 a and realizing that the potential of the ion 
with charge Zie is a static Coulomb potential in its rest 

frame, we find that both S _^ and _are at least of order 

Z 2 a. Using the recursion relation of Eq. (^^, the second 
term on the right hand side is at least of order {Z 2 a)^ and 
can therefore be neglected. Therefore, if only contributions 
at first order in Z 2 a are kept (whereas contributions at 
higher orders in Zia are important, as it is illustrated in 
Fig. 2), the matrix element for pair production is found 

to be A4 « 5-1_ S^_. This can also be seen directly from 

the definition A4 = _ SZ- and the fact that for a static 

field we have no off-diagonal elements in S with respect to 
the block structure given in Eq. ( |l5|) and therefore = 
S^_. Summing over all final states we get for Pi (keeping 
in mind that at this order of Z 2 a also \C\’^ = 1) 

Pi « Pi{S+-S+_S—S++) = Tr{S+-S^+). (49) 

Following the arguments for the high energy limit of Q| 
this probability should now be identical to the lowest or¬ 
der result. On the other hand the limit performed here 
corresponds exactly to the one which is done in order to 
derive the Bethe-Maximon theory ^,|^. Therefore the 
discrepancy between the two results is also present in this 
case, independently of the many-particle aspects. Com¬ 
paring the two multiplicities N and Npoisson one finds 
that they start to disagree only at fourth order in Z 2 a. 

Therefore one has to conclude that the discrepancies in 
the literature concerning Coulomb corrections cannot be 
traced back to the many-body aspects of pair production 
in this case. 
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Z, 



Z, 



Fig. 2. Examples of graphs contributing to pair production 
for Z 2 Q. <C 1. 

4 Conclusion and Outlook 

In this paper we have derived in a very straightforward 
way the link between the single-particle S matrix ele¬ 
ments obtained from the solution of the Dirac equation 
and its relation to the multiplicity and single-pair produc¬ 
tion probability. The term , which is missing in many 
papers and which is important for the calculation of the 
differential cross section with respect to both electrons 
and positrons, was introduced. A careful analysis of the 
limit which is relevant for a comparison with the Bethe- 
Maximon results shows that the absence of Coulomb cor¬ 
rections in this high energy limit cannot be understood by 
a purely conceptional treatment of the external field prob¬ 
lem. The problem therefore seems to lie in the way this 
high energy limit is made. That a too simplified use of the 
eikonal approximation leads to a result in contradiction 
with the Bethe-Maximon one was already pointed out in 
1 ^1, whereas within a more refined analysis as given in 
1^1 there is no disagreement. 

In priniciple the pair production in external fields is 
well understood, but the numerical evaluation of the Cou¬ 
lomb corrections for two highly charged ions, especially at 
small impact parameter, still needs to be found. 
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